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POTENTIAL FLOW SOLUTION FOR A STOL WING PROPULSION SYSTEM 


by James A. Albers and Merle C. Potter* 

Lewis Research Center 

SUMMARY 

The analysis considers a two-dimensional combined wing and propulsion system 
which consists of an airfoil (with flap) and fans located under the wing. 

A numerical method is developed which includes the effect of suction at the inlet of 
the propulsion system and treats the exhaust jet as part of the solid body to determine its 
approximate location. The method provides the potential flow solution for any fan or en- 
gine mass flow coefficient, the thrust coefficient for the propulsion system exhaust, and 
the wing and flap angles of attack. Validity of the numerical solution for a case with suc- 
tion (but with no jet) was indicated by application of the program to a two-dimensional in- 
let; excellent agreement was found with experimental results. 

The analysis was used to obtain pressure distributions, flow fields, and lift coeffi- 
cients for an externally blown flap, high-lift configuration. The flow field for this config- 
uration indicated high upwash angles (60° to 90°) at the propulsion system inlet and large 
jet penetrations at high angles of attack. A comparison of two-dimensional lift coeffi- 
cients obtained by the method of this report and Spence's jet flap theory indicated that the 
method of this report yielded lift coefficients that were an average of 10.5 and 12. 1 per- 
cent higher in the 30° and 60° flap angles, respectively. A comparison of three- 
dimensional lift coefficients with experimental data for the blown flap indicated good 
agreement for the 30° flap, with the predicted lift coefficient an average of 11.4 percent 
higher than experimental data. Calculated pressure distributions showed severe adverse 
pressure gradients over a large portion of the wing at angles of attack of 20° or greater. 


INTRODUCTION 

In recent years there has been much interest in short-takeoff -and-landing (STOL) 
aircraft for both civil and military applications. A STOL airplane must have the capa- 
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bility for both high lift at takeoff and low drag at cruise . Past experimental work (refs . 

1 to 3) demonstrated that the jet flap concept was capable of producing high lift. The jet 
flap airfoil injects high-velocity air over the flap surface from a slot located at the trail- 
ing edge of the airfoil, as shown in figure 1(a). (The jet flap airfoil (fig. 1(a)) is some- 
times referred to as a ”jet augmented flap” in the literature.) One way to implement the 
jet flap concept is to use an externally blown flap. This may be accomplished by using 
high -bypass -ratio turbofan engines which exhaust into the wing flap system. 

A STOL concept under investigation at Lewis is a multiple-fan externally blown flap 
(fig. 1(b)) . Important in this design concept is passing some of the fan exhaust through 
the gaps in the flaps to control the boundary layer over the wing upper surface. Some of 
the aerodynamic problems associated with this concept are (1) airfoil design for takeoff, 
cruise, and landing; (2) fan location and orientation; (3) jet penetration of propulsion 
system exhaust; and (4) the slot location and amount of blowing which are necessary for 
satisfactory boundary layer control. An analytical tool must be available to do detailed 
design studies of these aerodynamic problems. This tool must have the capability to 
handle both potential flows and boundary layers. 

The potential flow analysis is the first step in obtaining an analytical tool to design 
STOL wing propulsion systems. By shaping the airfoil geometry, the designer can mod- 
ify the wing pressure distribution to delay separation. Fan location and orientation can 
be improved by analysis of the pressure distributions and flow fields obtained from the 
potential flow solution of the combined wing and propulsion system. It is important to 
know the jet penetration of the propulsion system exhaust jet when considering the ver- 
tical location of the wing for a given aircraft configuration. Large penetration, under 
certain conditions, may cause a decrease in lift due to ground effect. Slot location and 
amount of blowing were not treated in this report since this would probably best be 
handled by boundary layer analysis (see ref. 4). From the potential flow solution, we can 
determine the maximum attainable lift coefficient for the wing propulsion system. Using 
the pressure distributions as inputs to a boundary layer program, the drag for a given 
engine-wing combination can be minimized. 

There are many approximate potential flow theories. Some approximate methods for 
calculating flow over two-dimensional bodies are discussed in references 5 to 7. Most 
approximate methods assume, for simplification, that the body is slender or that the per- 
turbation velocities caused by the body are small. Another type of approximate solution 
utilizes a distribution of singularities on or interior to the body surface. Some of these 
methods, based on a distribution of vorticity over the body surface are discussed in ref- 
erences 8 to 10. The potential flow theory that is often used when considering high-lift 
wing systems is Spence’s flap theory, as discussed in references 11 and 12. This thin 
airfoil theory considers the effect of the jet sheet leaving the trailing edge of the flap, but 
does not take into account the effect of the propulsion system inlet and the thickness dis- 
tribution of the lifting system. 
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The most general and comprehensive two-dimensional incompressible potential flow 
method and program is the Douglas method reported in references 13 to 15. This method 
utilizes a distribution of sources and sinks on the body surfaces, and does not require 
bodies to be slender nor perturbation velocities to be small. The program has the poten- 
tial for dealing with distributed suction over part of the surface, and hence, can handle 
the propulsion system inlet airflow. However, the program cannot handle problems for 
which the location of part of the boundary is unknown. For a combined wing and propul- 
sion system, the shape and location of the jet exhaust of the fan or engine is not known a 
priori; hence, a method is necessary to determine them. 

This report presents an analysis to solve the potential flow for a STOL wing propul- 
sion system, and the application of this analysis to a particular, externally blown flap, 
high-lift configuration. The analysis was accomplished by extending the two-dimensional 
Douglas method (ref. 14) to include the effect of suction at the propulsion system inlet, 
and by the development of a technique for determining the approximate location of the ex- 
haust jet of the propulsion system. This analysis was used to obtain pressure distribu- 
tions, flow fields, and lift coefficients for the externally blown flap configuration. 


ANALYSIS 

Representation of Wing Propulsion System 

While the present development can be used for any two-dimensional configuration, it 
is helpful in describing the analysis to consider a particular physical system. The high- 
lift wing propulsion system for STOL applications under investigation at Lewis is a 
multiple-fan externally blown flap, as shown in figure 2(a). The wind tunnel model is 
semispan with a NASA 4415 airfoil section, a 66-centimeter (26-in.) chord and a 165. 1- 
centimeter (65-in.) span. The model has eight propulsion units spaced spanwise with the 
inlets under the wing and the exhausts ahead of a double slotted flap. The 30° and 60° 
flap deflections in figure 1(b) represent typical takeoff and landing configurations. 

Since the proposed STOL lifting system utilizes a large number of fans closely spaced 
spanwise on each wing, it is reasonable to approximate the wing section flow character- 
istics with a two-dimensional flow. This approximation should be valid as long as there 
is a sufficient number of fans for blowing to be uniformly distributed along the wing 
trailing edge. The representation of the two-dimensional lifting system is shown in fig- 
ure 2(b). The equivalent body surface over which the potential flow is calculated con- 
sists of the airfoil with flap; the fans, which have a distributed suction at their inlet and 
a jet at their exit; and the jet sheet leaving the flap trailing edge. 

The wing - propulsion system combination is idealized by considering it to be one 
solid body with suction at the fan inlet. The jet stream, as it exists from the propulsion 
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system is at a higher total pressure than the surrounding flow. In potential flow the total 
pressure is everywhere constant; hence, in this study the jet is considered to be part of 
the solid body. This assumes no mixing of the external free stream and the free jet. The 
equivalent two-dimensional-propulsion-system dimensions and jet-sheet thickness were 
determined from the known mass flow rate and thrust of the I ewis propulsion system 
(fig. 2(a)). The method used to determine the location of the free jet is discussed in a 
later section. 

The potential flow problem for a given wing - propulsion system combination be- 
comes one of calculating the velocities on and external to the body surface for any com- 
bination of the following variables: (1) free-stream velocity V M , (2) fan or engine mass 
flow rate m per unit span, (3) propulsion system thrust T per unit span, (4) flap angle 
0, and (5) wing angle of attack a. The first three variables can be combined into two 
dimensionless parameters: the fan or engine mass flow coefficient Cq = m/pV^C and 

the thrust coefficient C T = t /( i / 2 pV 2 cj . The development of the theory to handle this 

calculation is discussed in the following sections. (All symbols are defined in appen- 
dix A.) 


Basic Equations and Boundary Conditions 

The basic potential flow equation is obtained from the incompressible continuity equa- 
tion together with the condition of irrotationality which gives LaPlace's equation 

V 2 <? = 0 (1) 

where cp is the velocity potential due to the presence of the body only. To ensure 
uniqueness of the solution, the regularity condition at infinity is specified as 

Moo- 0 (2) 

The velocity field V can be expressed as the sum of the two velocities 

V = + v (3) 

where is the free-stream velocity and v is the disturbance velocity due to the 
presence of the body only. 

A general method of solving the potential flow for an arbitrary boundary is by using 
a large number of sources and sinks distributed on the surface of the body. This is the 
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method presented in this report. The boundary condition, illustrated in figure 3, speci- 
fies that the entire normal component of velocity of the fluid at point p must be equal to 
the prescribed normal velocity on the surface. The contribution supplied by the source- 
sink distribution is v-n and that supplied by the free -stream velocity is V^-n. The 
prescribed normal velocity on the surface is due to suction or blowing. Thus the 
boundary condition becomes 


Voo-nlp -v.n| p = V N (4) 

Since v-n= dcp/dn, the boundary condition on cp is 

— - V N < 5 > 

an p 'p 

Equations (1), (2), and (5) form the classic Neumann problem of potential theory, which 
is the basic problem we wish to solve. The direct problem, as just defined, can be 
solved analytically by conformal transformation only for a limited class of boundary sur- 
faces. By using a large number of sources and sinks distributed on the surface of the 
body, the boundary condition can be formulated into an integral equation. 


Formulation of the Boundary Condition as an Integral Equation 

A simple potential function which satisfied equation (1) is the potential due to a point 
source. The potential at a point P due to source at q is expressed as 


d<P (P) = °<q> ds . 

r(P,q) 


( 6 ) 


where cr(q) is the local intensity per unit area of the source and r(P,q) is the distance 
between P and q. Because LaPlace T s equation is linear, the combined potential due to 
a distribution of sources is also a solution. By considering a continuous source distribu- 
tion on the surface S, the potential at point P due to the entire body becomes 


cp( P) = 



- ct M- ds 

r(P,q) 


(7) 
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The potential as thus given satisfies equations (1) and (2), but it must also satisfy the 
boundary condition as given by equation (5) . Applying the boundary condition requires 
evaluating the derivative dcp/d n at point p on the boundary surface . The derivative of 
l/r(p,q) becomes singular at p when p and q coincide, so that the principal value of 
the integral must be extracted. The principal value, according to reference 16, is 
-2rrcr(p). This is the contribution to the normal velocity at p from the source at p. 

The contribution of the remainder of the sources to the normal velocity is given by the 
derivative of the integral of equation (7) evaluated on the boundary. The normal deriva- 
tive of <p becomes 



( 8 ) 


Applying the condition of equation (5) to equation (8) results in the integral equation for 
the source -intensity distribution cr(p) 


2ircr(p) - 



(9) 


This equation is a Fredholm integral equation of the second kind whose solution is the 
central problem of the analysis. 

The quantity -3/3n j^l/r(p,q)j is called the kernel of the integral equation and de- 
pends only on the geometry of the surface. The first term of equation (9) is the normal 
velocity induced at p by a source at p. The second term is the combined effect of the 
sources at other points q on the surface of the body. The specific boundary conditions 
determine the right side of equation (9) . The first term on the right is the normal com- 
ponent of the free-stream velocity at p. The second term on the right is the prescribed 
normal velocity on the boundary surface at p. The solution of this Fredholm integral 
equation then requires determining the unknown function a on the body surface . 


Solution of Integral Equation 

Since the boundary of the wing propulsion system is completely arbitrary, the inte- 
gration of equation (9) with respect to S should be done numerically. The boundary is 
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approximated by a large number of surface elements whose characteristic lengths are 
small compared to the body. It is assumed that the surface element is a flat segment, as 
shown in figure 4. As the number of elements increases, the assumed model approaches 
the shape of the body. The value of the source intensity is assumed to be constant over 
each surface element. By assuming this constant intensity over each element, the prob- 
lem becomes one of finding a finite number of values of ex , one for each of the surface 
elements. This gives a number of linear equations equal to the number of unknown val- 
ues of cr . On each element a control point (the midpoint of the element) is selected 
where equation (8) is required to hold. Rewriting equation (8) in summation form yields 


d( p 

= -2uo(p) + \ — 

1 

9n 

P A 3 " 

r(p,q)_ 


p^q 


cr(q) AS 


( 10 ) 


The right side of equation (10) now becomes a matrix consisting of the normal velocities 
induced by a source of intensity a at the control points of all elements. The normal 

■Hi 

velocity at the control point of the i Ln element due to all surface elements is denoted as 


Thus 



A ii°i + 




A 


ij 


d 

i 

dn 

r(p,q)_ 


AS 


(ID 


where i corresponds to p and j corresponds to q. 

The source densities of all the surface elements must be determined in such a way 
that the normal velocity condition is satisfied at all control points. This results in 


N 


^ A..CX. 

/ , J 
j=l 



• n i + V N,i 


( 12 ) 


This set of linear algebraic equations is an approximation to the integral equation (9) . 
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Both V^-n and the prescribed normal velocity V^, in general, vary over the body sur- 
face. The linear equations are solved by a procedure of successive orthogonalization, 
as discussed in reference 17. Once the linear equations have been solved, flow veloci- 
ties may be calculated for points on and off the body surface . The method just described 
is used to obtain the basic solutions of potential flow. 


Basic Solutions 


The superposition of any solutions to the integral equation (9) is also a solution since 
LaPlace's equation is linear. Hence, the flow about a body may be thought of as a linear 
combination of four basic flows illustrated in figure 5: 

(1) Uniform flow at zero angle of attack 

(2) Uniform flow at 90° angle of attack 

(3) Vortex flow 

(4) Flow due to suction or blowing 

Uniform flow solutions . - The uniform flow solutions are solutions due to free- 
stream velocity (rectilinear flow) past the body surface at 0° and 90°, respectively. For 
these basic solutions, the boundary condition of zero normal velocity on the surface must 
be satisfied. Then the prescribed velocity normal to the surface must be zero for the 
basic uniform flow solutions. From equation (5) the boundary condition becomes 


d(p 

0n 


= V co’ n 


(13) 


The solution for the body at any angle of attack may be obtained by a linear combination 
of the 0° and 90° uniform flow solutions. 

Vortex flow solution . - For a lifting body the circulation is obtained by placing a 
vortex at any convenient location within the body. The boundary condition of zero normal 
velocity on the surface still applies (eq. (5)) except that now V^ is replaced by the vor- 
tex velocity at any point. If V y represents the velocity at any point p on the body 
caused by the vortex, the boundary condition for the basic vortex solution becomes 


d(£_ 

Sn 


= V v -n 


(14) 


Suction flow solution . - The suction flow solution is obtained by specifying a pre- 
scribed normal velocity V N on the surface of the body with a zero free -stream velocity. 
From equation (5) the boundary condition for the basic suction velocity solution becomes 
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(15) 


dc£_ 

3n 


-V 


N 


For each basic solution, the velocities on the body surface and at prescribed loca- 
tions in the flow field may be obtained. From the basic solutions the total combined 
solution may be obtained. 


Combination Solution 

The total velocity tangent to the body surface can be obtained by adding the tangential 
velocities of the four basic solutions. 

' V t = V t,0 cos a + V t,90 sin + rV t,v + V t, s (16 > 

where a is the angle of attack. 

The circulation T is determined by satisfying the Kutta condition at the trailing edge 
of the body. This condition stipulates that the flow at the body trailing edge be smooth. 
Thus, the tangential velocities above and below the trailing edge must be equal in magni- 
tude. If AV is defined as AV = V U pp er - V^ QWer , the Kutta condition is satisfied if 
AV = 0 at the body trailing edge. Then from equation (16) 

AV t>0 cos « + AV t 90 sin a + AV t s ♦ r AV tiV = 0 (17) 


Solving for T yields 


AV, n cos a + AV + Qn sin a + AV + c 

r = -- (is) 

A Vi. 

V t,v 

Once the combined velocities on the body surface are known, the pressure coeffi- 
cient, the lift coefficient, and the thrust coefficient can be obtained (see appendix B). 

For off -body points it is more convenient to combine the basic source intensities 
rather than the basic velocities. The equation for the combined source intensity is 

a = cr Q cos a + Og 0 sin a + ra v + a g (19) 

Then, the x and y components of velocity are calculated from the combined source in- 
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tensities. This approach is the same as that used in reference 15, with the addition of 
the basic suction source intensity being added to the other basic source intensities. 


Location of Propulsion System Exhaust Jet 

The location of the propulsion system exhaust jet is determined by the following var- 
iables: (1) jet angle at flap trailing edge 9 , (2) jet penetration H, (3) jet angle at trail- 
ing edge of jet 9 ^ , and (4) total length of the free jet L T . The representation of these 
variables is shown in figure 2(b) . It was assumed that the free jet leaves the flap trail- 
ing edge at the flap angle 9 . The flap angle is defined as the angle between the free- 
stream direction and lower flap surface. After the jet leaves the trailing edge, the free- 
stream velocity turns the jet, which approaches a horizontal asymptote at several chord 
lengths beyond the airfoil leading edge. For typical thrust coefficients Crp correspond- 
ing to takeoff and landing conditions, this occurs at approximately four chord lengths 
from the airfoil leading edge (see ref. 10). 

For a reasonable approximation to the jet shape, the lift coefficient is expected to 
depend principally on the vertical location of the jet asymptote. The problem then be- 
comes one of finding the jet penetration H as shown in figure 2(b). Initially, the pene- 
tration was assumed and a cubic equation used to approximate the shape of the jet sheet. 
A cubic equation is the simplest expression that adequately approximates the jet shapes 
obtained from Spence's jet flap theory of reference 12. The correct distance H is that 
value for which the vertical component of thrust at the flap trailing edge balances the in- 
tegrated vertical pressure forces on the free jet. Several values of H were assumed 
until the correct value of H was obtained. 

Since the angle of the free jet is not exactly horizontal several chord lengths beyond 
the wing, a small angle 9 ^ (5° or less) was assumed. The length of the jet Lrp was ex- 
tended until the vertical component of force on the end of the jet (last 5 percent of the jet) 
was negligible for the chosen angle 0j. Thus, if the jet is extended beyond this length, 
it gives no significant contribution to the lift coefficient. Neglecting the vertical compo- 
nent of thrust at the end of the jet for an angle of 5° results in only a 3 -percent variation 
in lift coefficient. Since the body representing the wing, flap, and jet is considered a 
lifting body, it was assumed to be closed at the jet trailing edge and the Kutta condition 
was applied here. 


Inlet Air Suction 

The inlet air suction was obtained by prescribing the velocity profile at the fan face. 
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From the prescribed normal velocity V N on the surface of the body, the basic suction 
solution is obtained. This gives the desired mass flow rate for the inlet of the propulsion 
system. 


DISCUSSION 

Validity of Analysis 

Inlet air suction . - To help ensure validity of the analysis, comparisons were made 
with known existing flow solutions. One existing solution solves the suction problem in- 
directly. It is also based on the Douglas method, but has application only to inlets and 
ducts (see ref. 18). This method utilizes three basic solutions, shown in figure 6, to ob- 
tain a combined solution of physical interest. The flow about the inlet is obtained by con- 
sidering the three basic solutions: Vj with inlet duct extension closed, Vg with the duct 
open, and the crossflow solution. With these three solutions any combination of 
free-stream velocity and mass flow through the inlet can be obtained. The duct must be 
extended far downstream of the region of interest to obtain valid solutions. This method 
could not be used to get solutions for a wing-engine combination since the body must be 
closed to consider it a lifting body. To make a comparison between this existing flow 
solution and the method presented in this report a two-dimensional inlet, shown in fig- 
ure 7, was considered. This inlet was chosen because experimental data were available. 
In the present analysis, mass flow through the inlet was obtained by considering a dis- 
tributed suction V N downstream of the inlet (see fig. 7). Comparison of the nondimen - 
sional surface velocities for the two methods are shown in figures 8(a) and (b) . Also 
shown is experimental data obtained from reference 19. The reference velocity V ^ 
was arbitrarily selected as the average velocity at an x/L of 0.89. Agreement between 
the two predictions is excellent for both the surface and centerline velocities. Compari- 
son of experimental data with the prediction is quite acceptable for the centerline veloci- 
ties. There is a slight variation between the experimental and predicted surface veloci- 
ties. One reason for this variation could be boundary layer effects. The preceding dis- 
cussion indicates that the combined uniform flow and suction solution is valid. 

Exhaust jet shape . - For a valid solution of a wing propulsion system there must be 
a reasonable approximation to the jet shape. The lift coefficients and pressure distribu- 
tions for a given thrust coefficient depend principally on the flap angle 6 and jet pene- 
tration H, as outlined in the section Location of Propulsion System Exhaust Jet, and not 
on the precise local shape of the jet. This is illustrated in figure 9, which considers 
various free jet shapes for a 30° flap. For clarity the jet thickness is not shown. The 
assumed cubic equation is shown, along with representative upper and lower bounds for 
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the jet shape. For the jet shapes A and C shown, the solution results in only a ±2. 5 per- 
cent variation in lift coefficient from the assumed cubic shape B. This percent variation 
is distributed over the entire wing surface, as illustrated by the pressure distributions 
in figure 9(b). Figure 9(b) shows only a 3 -percent variation in pressure distribution for 
the jet shapes considered. Thus, the lift coefficients and pressure distributions depend 
principally on the flap angle and jet penetration and not on the precise local shape of the 
jet for the present configuration. 

As a point of interest, a comparison of the jet shape based on the Spence's theory of 
reference 12 was made with the jet shape obtained from the present method. Refer- 
ence 12 assumes that all flow deflections from the free stream are not large and vortex 
distributions are placed on the x-axis rather than on the airfoil or jet. Thus, a compari- 
son could only be made for relatively small flap deflections (30° or less) . A comparison 
of the nondimensional jet shape predicted from Spence's theory and from the method of 
this report is shown in figure 10 for a 30° flap deflection and a thrust coefficient of 3. 

The basic shapes of the two cases are the same close to the wing. The jet penetration of 
the present method is larger than Spence's theory at the greater distances, as would be 
expected. The present method, besides not assuming small angle approximations, in- 
cludes the wing thickness and camber effect which would increase the lift coefficient and 
would also result in a greater penetration. 


Example Applications 

Flow field . - Potential flow solutions are adequate representations of the flow around 
bodies if the surface boundary layers are thin and remain attached. It is assumed that 
the final design of a high -lift wing propulsion system will be one in which boundary layer 
separation is prevented at relatively high angles of attack and flap settings. Representa- 
tive flow fields for an externally blown flap, high-lift configuration are shown in figures 
11(a) to (c). The flow fields were obtained by sketching streamlines tangent to the calcu- 
lated velocity vectors at various points in the flow. The wing propulsion system is 
shown, along with the shape of the jet exhaust of the propulsion system. For the condi- 
tions shown, the upwash angles at the propulsion system inlet are quite large, varying 
from 60° to 90° depending on flap angles and wing angles of attack. Two stagnation points 
occur on the lifting body. One occurs ahead of the inlet below the leading edge of the 
wing, and the other occurs downstream of the inlet on the under surface of the fan. Both 
stagnation points move further aft as the flap angle and the wing angle are increased. By 
observation of the flow fields it is seen that the under surface of the wing is in a rela- 
tively stagnant region. The jet penetration increases with angle of attack (figs. 11(a) 
and (b)). For a flap angle of 60° (fig. 11(c)) the jet penetration distance is approximately 
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three chord lengths at five chord lengths beyond the wing leading edge. This jet pene- 
tration is important when considering the effect of the ground on lift coefficient. 

Pressure distribution . - The predicted pressures on the surface of the airfoil are 
valid only if the boundary layer is very thin and attached to the surface. The potential 
flow pressure distributions can be used both to calculate the boundary layer growth on 
the surface of the airfoil and as a design aid for the combined wing and propulsion sys- 
tem. Pressure distributions on the wing upper surface with a 30° blown flap at various 
angles of attack are presented in figure 12. The incompressible pressure coefficient, 
corresponding to the minimum pressure point, ranges from -7. 5 to -51 for the 0° and 
20° angle of attack, respectively. For the high negative pressure coefficients, compres- 
sibility effects should be taken into account. The compressible pressure distributions 
can be obtained by making standard compressibility corrections to the incompressible 
pressure distributions (ref. 20). These extremely high negative pressure coefficients 
correspond to the very high lift coefficients which are discussed in the following section. 
The minimum pressure point for all angles of attack occurs very near the leading edge of 
the airfoil, and severe adverse pressure gradients over a large portion of the wing result 
at the higher angles of attack. The stagnation point moves further under the leading edge 
as angle of attack increases, resulting in high velocity gradients about the leading edge. 

To illustrate the effect of the inlet airflow of the propulsion system and the effect of 
the exhaust jet a comparison was made of the pressure distributions for (1) the wing alone, 
(2) the wing with jet but without inlet air suction, and (3) the wing with inlet air suction and 
jet. This comparison is presented in figure 13 for a 30° flap. At the minimum pressure 
point for the wing alone there exists a pressure coefficient of -4. 8 near the wing leading 
edge, followed by a mild adverse pressure gradient. The wing with jet but without inlet 
air suction would be representative of a jet flap airfoil shown in figure 1(a). Jet flap 
theory does not include the effect of the inlet airflow of the propulsion system. For the 
wing with jet (without suction) the pressure coefficient is about -18 at the minimum pres- 
sure point, and there is a severe adverse pressure gradient over a large portion of the 
wing upper surface. When the effect of the suction at the propulsion system inlet is in- 
cluded, the magnitude of the pressures is reduced considerably over the wing upper sur- 
face, resulting in a minimum pressure coefficient of -7. 6, followed by a much milder ad- 
verse pressure gradient. It may appear from the upper-surface pressure distributions of 
figure 13that theliftwith the jet alone is much larger than the lift associated with the jet 
with suction; but this is not the case if both upper and lower surfaces of the airfoil are 
considered. The change in pressure distribution between the zero suction case and the 
suction case is a result of a shift in the stagnation point (C = 1.0) on the under surface 
of the wing. For the wing without suction, one stagnation point occurs just ahead of the 
inlet of the propulsion system. For the wing with suction, this stagnation point moves 
closer to the wing leading edge and another stagnation point occurs on the under surface 
of the fan (fig. 11(a)). The corresponding shift in the pressure distributions on both the 
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upper and lower surfaces presented in figure 13 results in less than 5 percent decrease 
in lift when the effect of inlet suction is included for the selected inlet location. The pre- 
ceding discussion indicates that the effect of suction resulting from a fan or inlet in- 
stalled under the wing affects the pressure distribution on the wing upper surface favor- 
ably, with only a small effect on total lift coefficient. 

Lift coefficients . - In order to further indicate the applicability of the present anal- 
ysis a comparison (fig. 14) was made between Spence's theory (ref. 12) and the method 
of this report for two-dimensional lift coefficients for the blown flap configuration (fig. 
1(b)). The lift coefficients predicted by the method of this report for the 30° flap range 
from 6.5 to 13, while those for the 60° flap range from 15 to 21. The lift coefficients 
predicted by the present method generally range from 9. 1 to 12 percent and from 10.6 to 
13.6 percent higher than Spence's theory for the 30° and 60° flap, respectively. This 
difference exists since Spence's theory does not take into account the effects of the thick- 
ness and camber of the wing. The suction effect decreases the lift by approximately 
5 percent, as discussed previously. The thickness and camber effect corresponds to ap- 
proximately 15 percent increase in lift coefficient. 

The two-dimensional lift coefficients were used to determine three-dimensional lift 
coefficients to compare with experimental data of a semispan blown flap model (fig. 2(a)). 
The three-dimensional lift coefficient is 


C L = fCj (20) 

where f is a function of aspect ratio and thrust coefficient (assuming an elliptical lift 
distribution) and was obtained from reference 21 as 

2Crp 
AR + — - 

f = — - (21) 

AR + 2+0.604(C t ) 1//2 + 0. 876C T 

Calculated three-dimensional lift coefficients along with experimental data obtained from 
the Lewis test program are presented in figure 15. The aspect ratio was 5 for the blown 
flap model. The lift coefficients predicted by the method of this report range from 4 to 
7. 5 and from 9 to 13 for the 30° and 60° flap, respectively. The lift coefficients pre- 
dicted by the method of this report for the 30° flap (fig. 15(a)) range from 10. 8 to 12 per- 
cent higher than the experimental data. There is good agreement between Spence's 
theory and experiment for the 30° flap case. The lift coefficients predicted by the meth- 
od of this report for the 60° flap (fig. 15(b)) range from 28. 6 to 28. 4 percent greater 
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than the experimental data, while those predicted by Spence's theory range from 13. 5 to 
18. 1 percent greater than the data. 

The lift coefficient calculated from potential theory is the maximum attainable lift 
coefficient for each configuration corresponding to complete boundary layer control and 
negligible viscous effects. Thus, it is expected that the predicted lift coefficients be 
greater than the experimental data. The large variation in lift coefficients between 
theory and experiment for the 60° flap configuration may indicate that this configuration 
did not have optimum boundary layer control and that improvements could be made in 
obtaining better experimental coefficients. 


CONCLUDING REMARKS 

A method was developed to determine the two-dimensional potential flow solution of 
STOL wing propulsion systems. The Douglas potential flow analysis was extended to in- 
clude the effect of suction at the propulsion system inlet and to provide a technique for 
determining the approximate location of the exhaust jet of the propulsion system. The 
effect of suction was obtained by combining the Douglas basic suction solution with the 
uniform flow solution for a lifting body. The jet exhaust was considered as part of the 
solid body and its location was determined by balancing the vertical component of thrust 
at the flap with the integrated vertical pressure forces on the free jet. 

The applicability of the program is illustrated by considering a multiple -fan exter- 
nally blown flap under high-lift conditions. The results indicated high upwash angles at 
the fan inlet and large jet penetration at high angles of attack. The predicted two- 
dimensional lift coefficients for a 30° flap were an average of 10.5 percent higher than 
predicted by Spence's jet flap theory which neglects thickness effects. The predicted 
three-dimensional lift coefficients were an average of 11.4 percent higher than experi- 
mental data for the 30° blown flap high-lift configuration. The calculated pressure dis- 
tributions indicated that the minimum pressure point is near the leading edge (less than 
2 percent of chord) of the airfoil, with severe adverse pressure gradients at high angles 
of attack. 

The ability of the potential flow solution to predict pressure distributions, lift coef- 
ficients, and flow fields makes it extremely useful as a tool in the analysis and design of 
a STOL -wing propulsion system. The pressure distributions obtained from the potential 
flow solution can be used to determine the boundary layer growth and separation on the 
airfoil of the wing propulsion system. These distributions can be used to design the air- 
foil and to determine the optimum location and orientation of the propulsion system. 

From the predicted pressure distributions and boundary layer calculations, the frictional 
drag can be obtained for a given wing propulsion system. The potential flow solution can 
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be used to determine the jet penetration - an important quantity when considering ground 
effect. 

The analysis has application not only to wing propulsion systems, but to any lifting 
or nonlifting body where suction or blowing is applied. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 7, 1971, 

721-03. 
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APPENDIX A 


SYMBOLS 

Ay normal velocity of element i caused by a unit source at element i 

A^. normal velocity of element i caused by a unit source at element j 

AR aspect ratio 

C chord length 

three-dimensional lift coefficient 
two-dimensional lift coefficient 
C pressure coefficient 

r 

Cq mass flow coefficient 

C,p thrust coefficient 

F force in vertical direction 

f correction factor (eq. (21)) 

H jet penetration (fig. 2(b)) 

L length 

L,p total length of free jet 

M number of elements that describe the jet 
m fan or engine mass flow rate per unit span 
N number of elements that describe the body 
n a normal to the body surface 

P ' arbitrary point in the flow field off the surface 

p static pressure 

r distance between two points 

S surface of body 

T thrust per unit span 

V velocity 

v disturbance velocity 

x Cartesian coordinate 

y Cartesian coordinate 
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a angle of attack 

orientation of surface element 
r nondimensional circulation 

7} turning efficiency of exhaust jet 

8 flap angle (fig. 2(b)) 

jet angle at trailing edge of jet (fig. 2(b)) 
a surface source intensity per unit area 
p density 

c p velocity potential 

Subscripts: 

j-i_ 

i control point of i element 

ai_ 

j control point of j element 

N normal 

p arbitrary point on the surface 

q a surface point 

ref reference 

s suction flow solution 

t tangential 

v vortex flow solution 

00 free stream 

0 flow solution at zero angle of attack 

90 flow solution at 90° angle of attack 
Superscript: 

— vector 
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APPENDIX B 


CALCULATION OF FLOW QUANTITIES 

Once the combined velocities on the body surface are calculated, the pressure dis- 
tribution and the lift coefficient of the body can be found. The equation of motion for 
steady, incompressible, inviscid fluid can be expressed as 

(V. V)V = Vp (Bl) 

P 


For potential (ir rotational) flow, Bernoulli's equation results 


£ + i = Constant 
P 2 


(B2) 


The pressure coefficient Cp is defined as 


„ P - Poo 
Cp _ 1 2 


(B3) 


By use of equation (B2) 


c = i - 
P 2 

* oo 


(B4) 


The two-dimensional lift coefficient is defined as 


c z = 


2 


(B5) 


This can be obtained by integration of the pressure distribution over the surface of the 
body. 


r . Since L = / Pj cos a. dSj, 

•/s 
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where T is the exit thrust of the propulsion system. The exit thrust is obtained from 
the vertical force on the jet, the jet deflection angle, and the experimental turning effi- 
ciency rj between the propulsion system exhaust and the trailing edge of the flap. Then 


F 

T = l (B8) 

(sin 6)(ij) 

where r\ is the turning efficiency of the exhaust jet. The vertical force is calculated by 
integration of the pressures on the jet 

M 

F y = P i COS AS. (B9) 

i=l 

where is the angular orientation of the 1 element. 
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(b) Two-dimensional representation of wing propulsion system 
Figure 2. - Concluded. 


P. 



Figure 3. - Representation of boundary condition on 
body surface. 
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Jet shape 



(b) Upper-surface pressure distributions. 

Figure 9. - Effect of jet shapes on upper-surface pressure distribution. Flap angle, 30 P; wing angle of attack, 0°. 



Figure 10. - Comparison of theoretical nondimensional jet shapes. Flap angle, 30 P; thrust 
coefficient, 3. 
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(b) Flap angle, 60 P. 

Figure 14. - Comparison of theoretical two-dimensional 
lift coefficients for blown flap. Thrust coefficient, 3. 
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(b) Flap angle, 60°. 

Figure 15. - Comparison of calculated and experimental 
three-dimensional lift coefficients for blown flap. 
Thrust coefficient, 3. 
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